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ABSTRACT 

We present the first interferometric polarization map of the W51 el/e2 molecular cores obtained with 
the BIMA array at 1.3 mm wavelength with approximately 3" resolution. The polarization angle varies 
smoothly across the double cores with an average position angle of 23°±5° for W51 el and 15°±7° for 
W51 e2. The inferred magnetic field direction is parallel to the minor axis of the double cores, which is 
consistent with the theoretical picture that clouds collapse along the field lines. However, the magnetic 
field may not determine the axis of angular momentum of these two cores as the field directions of the 
two cores significantly differ with the previously measured directions of rotational axes. The polarization 
percentage decreases toward regions with high intensity, suggesting that the dust alignment efficiency 
decreases toward high density regions. The field directions are highly ordered, and the small dispersion 
of the polarization angles implies that magnetic fields are strong (> 1 mG) and perhaps dominate 
turbulence in W51 el/e2. 



Subject headings: ISM: magnetic fields - ISM: individual: W51 
techniques: interferometric 



polarization - star: formation 



1. INTRODUCTION 

Star formation in molecular clouds is one of the most 
fundamental problems in astrophysics. It has become in- 
creasingly evident that star formation cannot be fully un- 
derstood without considering magnetic fields (see recent 
reviews by Mouschovias & Ciolek 1999, and Shu et al. 
1999). Magnetic fields provide explanations for the sup- 
port of clouds against self-gravity, the formation and evo- 
lution of cloud cores, the origin of supersonic line widths, 
the low specific angular momentum of cloud cores and 
stars, and the formation of bipolar outflows. Unfortu- 
nately the magnetic field remains the most poorly mea- 
sured quantity in the star formation process. In order to 
advance our understanding of star formation, it is essential 
to improve our empirical knowledge of magnetic fields in 
molecular clouds. 

The most promising technique to probe magnetic field 
morphology in dense molecular cores is to measure the 
linear polarization of the thermal emission from spinning, 
magnetically aligned dust grains (Heiles et al. 1993). In- 
terstellar dust grains are elongated with their minor axes 
aligned with the magnetic field, thus the magnetic field 
direction in the plane of sky is perpendicular to the di- 
rection of polarization (Davis & Greenstein 1951; Roberge 
1996). Millimeter interferometers detect the dust emission 
in the low optical depth regime with high angular resolu- 
tion, thereby revealing the field morphology in the dense 
cores. In the past few years, linear polarization observa- 
tions with the Berkeley-Illinois-Maryland Array (BIMA) 
of Orion-KL and NGC 1333 IRAS 4A have demonstrated 
the capability and reliability of polarization observations 
with this telescope (Rao et al. 1998; Girart et al. 1999). 
These results are not only consistent with the previous 



single-dish observations, but also provide new information 
at resolutions up to 3". Therefore, it is reasonable for us to 
carry out a large linear polarization survey of star-forming 
cores with BIMA. 

W51 is a large molecular cloud/Hn region complex at 
a distance of 7. Oil. 5 kpc (Genzel et al. 1981). Associated 
with ultracompact (UC) Hn regions (Scott 1978), strong 
infrared and submillimctcr continuum emission (Thron- 
son & Harper 1979; Jaffe, Becklin, & Hildebrand 1984), 
molecular gas (Ho, Genzel, & Das 1983) and H2O and OH 
masers (Genzel et al. 1981; Gaume & Mutel 1987), W51 
is one of the most active high-mass star formation regions 
in our Galaxy. W51 el and W51 e2 are molecular cores 
located in the eastern edge of W51. The mass of W51 el 
and e2 are ~ 150 and 110 Mq, derived from NH3 obser- 
vations (Zhang & Ho 1997). With infall motions identi- 
fied, these molecular cores are in the process of forming 
OB associations (Ho & Young 1996). Subarcsecond VLA 
observations reveal that W51 el contains three UC H11 
regions, el, e3, e4 (Gaume, Johnston, & Wilson 1993), 
and a possibly dust source, e8 (Zhang & Ho 1997). In the 
vicinity of the UC Hn region e2, there is a compact water 
maser concentration known as W51 Main (Genzel et al. 
1981). Previous single-dish measurements with resolution 
~30" show very low polarization in W51 el/e2 consistent 
with non-detection (Kane et al. 1993; Dotson at al. 2000); 
however, this could be caused by averaging over the large 
beams. Our observations with 10 times better resolution 
will allow us to explore the field structure in the W51 el/e2 
cores. 

In this paper, we present the first interferometric po- 
larization map of the W51 el/e2 cores at 1.3 millimeter 
wavelength. Other sources in our survey will be reported 
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in later papers. 

2. OBSERVATIONS AND DATA REDUCTION 

Several observations were carried out from 1999 August 
to 2000 April using nine BIMA antennas with 1-mm SIS 
receivers and quarter-wave plates. The digital correlator 
was setup to observe the continuum emission with a 750 
MHz window centered at 226.9 GHz in the lower sideband 
and a 700 MHz window centered at 230.9 GHz in the up- 
per sideband. Strong CO J=2-l line emission was isolated 
in an additional 50 MHz window in the upper sideband. 
The primary beam is ~ 50" at 1.3 mm wavelength. Data 
were obtained in the B, C and D array configurations with 
projected baseline ranges from 6-170 kilowavelengths, re- 
sulting in a linear polarization map of the W51 el/e2 cores 
with a long integration time, ~ 25 hours. 

The BIMA polarimeter and the calibration procedure 
are described in detail by Rao et al. (in preparation). A 
pair of quarter-wave plates was placed in front of the lin- 
early polarized receiver of each antenna in order to enable 
the detection of the left (L) and right (R) circularly polar- 
ized radiation. Four cross correlations, LL, RR, LR, and 
RL, must be measured at each baseline in order to obtain 
the four Stokes parameters, /, Q, U, and V (Thompson, 
Moran, & Swenson 1986). Because there was only one 
receiver at each BIMA antenna, the quarter-wave plates 
were switched on a time scale shorter than the uv-cell tran- 
sit time to achieved quasi-simultaneous dual polarization 
observations. Walsh function switching patterns (Thomp- 
son, Moran, & Swenson 1986; Harmuth 1969) were used 
to maximize the efficiency of obtaining four cross correla- 
tions (Rao 1999). Sixteen switch patterns were needed to 
complete a cycle. The integration time of each pattern was 
set to 11.5 seconds and the plate switch took an additional 
2-3 seconds; therefore, one cycle was about 4 minutes. 

Calibration and data reduction were carried out using 
the MIRIAD software (Sault, Tcuben, & Wright 1995). 
The instrumental polarization response ("leakage") was 
calibrated by observing the strongly polarized quasars 
3C273 or 3C279 over a wide hour angle range, typically 
more than 5 hours. The observed polarization, including 
contributions from the instrumental polarization and the 
polarization of the calibrator, varied with hour angle due 
to the fact that the polarization vector of the calibrator 
rotates with respect to the linear feed horn. On the other 
hand, the instrumental polarization is constant with time. 
Therefore, the leakages can be solved by assuming that the 
polarization of the calibrator was constant in the frame of 
sky over the whole track. Task GPCAL of MIRIAD was 
used to solve simultaneously for the leakages and antenna 
gains from the polarization calibrator. The average leak- 
age of each antenna was 5.9%±0.4% for our observations. 
The leakage correction of each antenna was applied to vis- 
ibility data before further processing. 

Channels with line emission were carefully flagged out in 
the visibility data of the continuum bands. As the contin- 
uum emission of the the W51 el/e2 cores was stronger than 
the phase calibrator (QSO 1751+096), self-calibration was 
performed for refining antenna gain solutions. The Stokes 
I image was made with Briggs' robust weighting of 0.5 
(Briggs 1995; Sault & Killeen 1998) to acquire a smaller 
beam size (2'.'7x2'.'0, PA=1°) for better determination of 



the self-calibration model. The model was used to cali- 
brate the visibility data in order to obtain new gain so- 
lutions with shorter time intervals. The above iterations 
were repeated until the Stokes / image was not sufficiently 
improved. The Stokes Q and U images were then made 
with natural weighting to obtain the highest S/N ratio. 
The resulting synthesized beam was 3'.'2x2'.'3 with PA=2°. 
Maps of Stokes Q and U were deconvolved and binned to 
approximately half-bcamwidth per pixel (r75xl'/2) to re- 
duce oversampling in our statistics. We then combined the 
maps to obtain the observed linearly polarized intensity, 

i P ,ob S = VQ 2 + u 2 , (l) 

and the polarization position angle, 

<j) = 0.5 tan" 1 (U/Q). (2) 

As the polarized intensity is an intrinsically positive quan- 
tity, Eq. (1) will tend to overestimate the polarized inten- 
sity. Therefore, a bias correction for the polarized intensity 
must be performed (Leahy 1989), 
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where the rms noise value of polarized intensity, aj , was 
taken as the average value of the rms noise in the Q and 
U maps, as the noise levels in these two maps are compa- 
rable. The polarization percentage was calculated from 



P= T- 



(4) 



The measurement uncertainty in the position angle, 

<j<p = 0.5 tan -1 (cr Ip /Ip), (5) 

depends on both oj and I p , thus it varies across the 
map. When weighted with I p , the average measurement 
uncertainty in the position angle for our observations was 
5?6±1?7. 

3. RESULTS AND ANALYSIS 

Figure 1 displays a contour map of the 1.3-mm contin- 
uum from the W51 el/e2 cores overlaid with polarization 
vectors. Polarization is called detected at positions where 
the linearly polarized intensity is greater than 3<7j (lcrj 
= 4.7 mjybeam -1 ) and the total intensity is greater than 
5cr/ (lcrj = 27 mjybeam -1 , which is dominated by in- 
complete deconvolution rather than thermal noise) . Under 
these criteria, the polarized emission extends over an area 
of ~12 beam sizes. Table 1 lists the polarization measure- 
ments in the W51 el/e2 cores at selected positions, sep- 
arated by approximately the synthesized beamsize. The 
distributions of the polarization angle and the polariza- 
tion percentage are shown in Figure 2 and 3, respectively. 

3.1. Continuum and Polarized Emission 

The south and north continuum condensations of the 
double cores in Figure 1 correspond to W51 el and W51 
e2, respectively. W51 el is extended (FWHM=5'.'0x3'.'0, 
PA=11°) and its integrated flux is 8.2 Jy. W51 e2 is also 
resolved (FWHM=4'/4x3'/3, PA=-8°) with an integrated 
flux of 13.6 Jy. We will refer to the intensity maxima of 
these two cores as "el 1mm peak" and "c2 1mm peak" . 
Analysis of the continuum flux at multiple wavelengths 
has shown that W51 el and W51 e2 contain both free-free 



emission and dust emission (Rudolph et al. 1990). The ex- 
pected flux of the free-free emission at 1.3 mm wavelength 
for these two cores can be estimated from their SEDs 
(Rudolph ct al. 1990), which are ~ 0.5 Jy for both cores. 
Therefore, the average fractions of the free-free emission 
in W51 el and e2 are ~6% and <~4%. 

The polarized emission does not show two distinct peaks 
associated with W51 el and e2, but it mainly arises from 
an extended region across the two cores. The peak of po- 
larized intensity is at l'/2 south of el 1mm peak (Table 
1). There is also a compact region of polarized emission 
^4" northwest of e2 1mm peak (for convenience, we name 
this position "e2 pol NW"; cf. Fig. 1). The polarized flux 
drops to zero between e2 1mm peak and e2 pol NW, and 
e2 pol NW shows very different polarization angle. There- 
fore, the polarization in e2 pol NW may sample a different 
region, and is excluded in our statistical analysis. 

Most of the polarized flux is associated with W51 el. 
The elongated shape of W51 el is due to the fact that it is 
composed of emission from nearby UC Hll regions (el, e3 
and e4), a newly discovered continuum source (e8: Zhang 
& Ho 1997), and the dust emission in the W51 el core. 
The nearest source to el 1mm peak is e8, which indicates 
that e8 is stronger than the UC Hll region el at 1.3 mm. 
Because the UC Hll region el is stronger than e8 at 1.3 
cm (Zhang & Ho 1997) and e8 is not detected at 2 cm 
(Gaume et al. 1993), the spectral index of e8 is very dif- 
ferent from free-free emission. This fact is consistent with 
Zhang & Ho's suggestion that e8 is a dust-dominated con- 
tinuum source. Hence, c8 is the most likely site for active 
star formation in the W51 el core. 

The position of the UC Hll region e2 (accurate to 0'.'2; 
Gaume ct al. 1993) is offset ~1" to northwest of e2 1mm 
peak . It is known that the molecular gas in W51 e2 is 
not evenly distributed around the H n region (Zhang & 
Ho 1997). If we consider that molecular gas is in gen- 
eral associated with dust grains, then, although the offset 
is small, it is evident that the free-free emission and the 
dust emission in W51 c2 do not coincide. Since the polar- 
ized intensity only originates from the dust emission and 
the fraction of the free- free emission could be significantly 
higher than the average value at the position of the UC 
H II region, this offset may be part of the reason for the 
absence of polarization between e2 1mm peak and e2 pol 
NW. It is interesting that the location of W51 Main is right 
in the polarization gap between e2 1mm peak and e2 pol 
NW, and a stream of water masers in W51 Main has been 
detected with proper motions headed toward position an- 
gle 200° (Leppanen, Liljestrom, & Diamond 1998), which 
seems well matched with the direction of the polarization 
gap. 

3.2. Polarization Angle Distribution 

The polarization angle distribution in the W51 el/e2 
cores appears to be well matched to the morphology of 
the cores. Except that e2 pol NW has a very different 
position angle at 58°±8°, the rest of the polarization vec- 
tors are uniformly distributed along the major axis of the 
double cores with an average position angle of 21°±6°. 
Polarization vectors seem to be approximately parallel to 
local core edges: in W51 el, polarization angles decrease 
westward and increase eastward from the central ridge to 



the edge; in W51 e2, position angles are close to 0° at 
the northeast edge; even for e2 pol NW, the difference be- 
tween the polarization angle and the edge is less than 20°. 
The polarization angles also follow the ridge connecting 
W51 el and W51 e2. Excluding e2 pol NW, we calcu- 
late the mean polarization angles for W51 el and W51 
e2 weighted by the measurement uncertainty, which are 
23°±5° for W51 el and 15°±7° for W51 e2. Compared 
to the measurement uncertainty, the variation of polariza- 
tion angles corresponding to the core morphology is only 
slightly suggested by our data. 

Figure 2 shows the distribution of the polarization an- 
gles in W51 el and W51 e2. The observed polarization 
angle dispersion S(f> b s , defined to be the standard devi- 
ation of the position angles, is 4?8±2?0 in W51 el and 
6?6±3?0 in W51 e2 (excluding e2 pol NW). Interestingly, 
these values are equal to the measurement uncertainty of 
the position angles a^, 4?9±1?4 in W51 el and 6?7±3?4 
in W51 c2. The observed dispersion is made up of contri- 
butions from the measurement uncertainty and the intrin- 
sic dispersion 5<j), S(f> 2 bs = 5(f) 2 + a?, which indicates that 
the intrinsic dispersions in W51 el and e2 are very small. 
Therefore, if we take o^ as small as possible (a^ — lo~o-,) 
and 5<f) bs as large as possible (6<p b s + l°~5<f> obs )i we can 
obtain the upper limit of the intrinsic dispersion, which 
is 5?8 in W51 el and 9?0 in W51 e2. Our results for the 
dispersion analysis are summarized in Table 2. 

3.3. Polarization Percentage Distribution 

The average polarization percentage of the W51 el/e2 
cores is 1.8%±0.4%, including positions for which polar- 
ization vectors are not plotted in Fig. 1 because I p is 
too small. W51 e2 has a lower polarization percentage 
(1.1%±0.3%) than that of W51 el (3.2%±0.6%). To bet- 
ter compare with previous observations, we convolved our 
map to 30" resolution with a Gaussian beam and ob- 
tained 1.3%±0.2% polarization with a position angle of 
25°±3°. The polarization detections we have obtained are 
clearly different from the single-dish results: at 1.3 mm 
(HPBW=30"), Kane et al. (1993) measured a polarization 
percentage of 0.50%±0.21% at an angle of -31°±11°; at 
100 ^m (HPBW=35"), Dotson at al. (2000) obtained zero 
polarization (0.54%±0.58%). However, all results could 
be still consistent, because the low S/N ratio of the these 
single-dish measurements can lead to large uncertainty in 
the polarization angle. Our observations demonstrate that 
high resolution is important for successful detection and 
mapping of magnetic fields in molecular cores. 

Figure 3 shows plots of the polarization percentage ver- 
sus Stokes /. These plots show that the polarization per- 
centage decreases as the total intensity increases in both 
W51 el and W51 e2. Such a decrease in the polarization 
toward high intensity regions (the depolarization) has been 
commonly seen in polarization observations, e.g., L1755 
(Lazarian, Goodman, & Myers 1997), OMC-1 (Schlcuning 
1998), and OMC-3 (Matthews & Wilson 2000). Lazar- 
ian, Goodman, & Myers (1997) have shown that the de- 
polarization can be explained as a consequence of the dust 
alignment efficiency decreasing toward the inner parts of 
the dark clouds, as all known mechanisms of grain align- 
ment fail under the typical physical conditions of dark 
cloud interiors. On the other hand, beam smearing over 



small-scale field structure can also produce low polariza- 
tion percentage (Rao et al. 1998); this can only be verified 
by observations with higher resolution. 

Due to the apparent close anti-correlation of the polar- 
ization percentage and the total intensity in Figure 3, it 
is interesting to perform a least-squares linear fit on these 
two quantities. We obtain log 10 (p) = (—1.52 ± 0.01) — 
(0.63 ± 0.07) x log 10 (7) with a correlation coefficient of 
-0.91inW51 el, and log 10 (p) = (-1.73 ± 0.01) - (0.99 ± 
0.04) x log 10 (7) with a correlation coefficient of —0.97 in 
W51 e2. However, due to the difficulty of excluding the 
free-free emission at every position, direct comparison be- 
tween our results and the dust alignment mechanisms can- 
not be done. The tight correlations hint that the depolar- 
ization in the W51 el/e2 cores is a gradual and global 
effect. Therefore, it is unlikely that the depolarization is 
entirely caused by a sudden change of the alignment mech- 
anism in a local region, as was reported in Orion-KL (Rao 
et al. 1998). 

4. DISCUSSION 
4.1. Magnetic Field Morphology 

The magnetic field direction inferred from the param- 
agnetic relaxation of grain alignment is perpendicular to 
the direction of the linear polarization of the dust emis- 
sion (Davis & Grccnstcin 1951). Therefore, the average 
directions of the magnetic fields are 113°±5° in W51 el 
and 105°±7° in W51 e2 (excluding c2 pol NW). The field 
directions are approximately parallel to the minor axis of 
the two cores, suggesting that matter collapsed along the 
field lines to form the parent cloud of the double cores. 
The uniform field structure in the W51 el/e2 cores also 
suggests that ionized gas in the UC Hn regions does not 
have much interaction with the magnetic field, or the inter- 
action only exists on a scale smaller than our synthesized 
beam. 

Comparison between the field morphology and the ro- 
tational axes of W51 el and W51 e2 provides interesting 
implications for the relation between magnetic fields and 
core formation. Rotational motions have been observed 
toward these two cores by Zhang, Ho, & Ohashi (1998) 
from analysis of the velocity gradient of CH3CN emission, 
and the derived rotational axis is at 66°±27° for el and at 
20°±17° for e2. Clearly, the rotational axes and the mag- 
netic field directions of these two cores do not align with 
each other. Although it has been suggested that protostel- 
lar disks may drag the field lines as they rotate (Holland et 
al. 1996), it may not be the case for these two cores. The 
field directions vary smoothly along the two cores. Thus, it 
is unlikely that two independent rotating cores that twist 
the field lines separately would produce a matched field 
morphology. Therefore, the smooth magnetic field struc- 
ture is most likely associated with the common envelope 
of el and e2. Zhang, Ho, & Ohashi (1998) caution that 
the uncertainty in their analysis may be large. However, 
if the difference between the directions of rotational axes 
and magnetic fields is indeed significant, the magnetic field 
of the parent cloud seems to have had no effect in deter- 
mining the axis of angular momentum of these two cores. 

4.2. Estimation of the Magnetic Field Strength 



Owing to the lack of understanding of the detailed mech- 
anism of grain alignment, the magnetic field strength can- 
not be directly derived from linear polarization measure- 
ments of dust emission (Lazarian, Goodman, & Myers 
1997). However, Chandrasckhar & Fermi (1953, hereafter 
CF) have proposed a method to estimate the field strength 
from the dispersion of polarization angles (<5</>) by assum- 
ing that the variation in polarization angles results from 
the perturbation of Alfvcn waves on the field lines. In this 
case, the field strength projected in the plane of the sky 
(B p ) is given by 

B p = Q^Tp %, (6) 



5<t> ' 

where p is the average density, 8vi os is the rms line-of-sight 
velocity, and Q is 1 for CF's derivation. Detailed studies 
show that the CF formula tends to overestimate the field 
strength. Q can be reduced by several factors, such as 
the inhomogencity of clouds and the line-of-sight averag- 
ing (Zweibel 1990; Myers & Goodman 1991). Ostriker, 
Stone, & Gammie (2001) calculate the value of Q from 
their simulations of turbulent clouds, and suggest that the 
CF formula with Q ~ 0.5 can account for the complex 
magnetic field and density structure, and provide accurate 
estimates of the plane-of-sky field strength under strong 
field cases when 5(j) < 25°. Since the angle dispersion we 
measured in the W51 el/e2 is small (5(f) < 9°; §3.2), it is 
interesting to work out the field strength estimates with 
the modified CF formula. 

However, the two other quantities, p and 8vi os , needed 
to complete this exercise have uncertainties not smaller 
than that of 5<fi. The density of the W51 el/e2 cores 
varies with the physical components and scales, and we 
should use the average density obtained from the region 
that the observed polarization is associated with. Zhang 
& Ho (1997) measure the NH3 emission in the infall re- 
gions of the cores and derive nn 2 ~ 3 x 10 6 cm~ 3 , which 
is probably too high for our purpose because the polarized 
intensity in the W51 el/e2 seems to originate in a more ex- 
tended region. From spectral decomposition of the whole 
W51 A region (including W51 el/e2 and W51 IRS 1), Sicv- 
ers et al. (1991) show that W51 A consists of a 20 K cold 
dust component with nn ~ 5x 10 4 cm~ 3 (n# = 2uh 2 ) 
and a 60 K warm dust component with nn ~ 10 6 cm~ 3 . 
This provides a range for the density as W51 el/e2 may 
contain a mixture of both components. Because the col- 
umn density of the warm dust is ~5 times more than that 
of the cold dust and W51el/e2 is the densest region in 
W51A (Sievers et al. 1991), W51el/e2 is more likely to 
have a density close to 10 6 cm~ 3 . It is also difficult to 
decompose Svi os that is associated with Alfvcnic motion 
from the molecular linewidths which are broadened by the 
complicated dynamics in W51 el/e2, such as infall (3.5 
km s _1 : Zhang et al. 1998), rotation (<~4 km s _1 : Zhang 
& Ho 1997), and possible outflow activities. The lincwidth 
of an optically thin NH3 line measured in the envelope of 
el and e2 by Young, Keto, & Ho (1998), Au=1.25 km s - \ 
may contain little contamination from dynamical motions. 
Therefore, we adopt Svi os — Au/v8hi2 = 0.53km s~ . 

Given the large uncertainty of the input parameters, we 
can only calculate the lower limits of B p using the upper 
limits of the angle dispersions obtained in §3.2. The de- 
rived B p for the density limits of W51 el/e2 are listed in 



Tabic 2. Our results show that the lowest value of B p in 
W51 el/e2 is just comparable to the typical magnetic field 
strengths measured in molecular cores with uh 2 = 10 5-6 
cm~ 3 using the Zeeman effect (<~ 500/iG: Crutcher 1999). 
Since this typical field strength provides crucial contribu- 
tion in cloud evolution (Crutcher 1999), our results show 
that the magnetic energy is an important component in 
the energetics of the W51 el/e2 envelope. It can be shown 
that the ratio of the turbulent to magnetic energy is simply 
proportional to the square of the angle dispersion (Lai et 
al. 2001, in preparation). The small angle dispersions we 
measured here imply that the magnetic field dominates the 
turbulent motion in the regions that emit polarized flux in 
the W51 el/e2 molecular cores. 

5. CONCLUSIONS 

We have measured linear polarization of the thermal 
dust emission at 1.3 mm toward the W51 el/e2 cores with 
the BIMA array. The conclusions of our observations are 
the following: 

• The magnetic field is parallel to the minor axis of 
the W51 el/e2 cores, consistent with the theoretical 



picture that cloud cores form due to gravitational 
contraction along field lines. However, the signifi- 
cant difference between the directions of magnetic 
field and rotational axes hints that the magnetic 
field of the parent cloud does not determine the 
axis of the angular momentum of these two cores. 

• The polarization percentage decreases toward the 
high intensity regions. The tight anti-correlation 
between log p and log I in both W51 el and W51 
e2 implies that depolarization is a global effect and 
may be caused by the decreasing dust alignment 
efficiency toward higher density regions. 

• The small dispersion of the polarization angles in 
W51 el/e2 suggests that the magnetic field is strong 
(> 1 mG) and dominates turbulence in the regions 
where the polarized flux arises. 
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Table 1 
Polarization Measurements in the W51 e1/e2 regions 



Position 


Stokes / 


Polarization 


Polarization 


Note 


III //\a 


( Jybeam -1 ) 


Percentage(%) 


Angle (°) 




(-6.3, 3.3) 


0.25±0.03 


6.4±2.0 


58± 8 


e2 pol NW 


(-1.0, 2.0) 


0.35±0.03 


5.6±1.4 


1± 7 


e2 


(-3.6, 0.6) 


3.22±0.03 


0.7±0.2 


19± 6 


e2 1mm peak 


(-1.6, 0.3) 


0.91±0.03 


1.5±0.5 


10±10 


e2 


(-1.6,-2.9) 


0.32±0.03 


7.9±1.6 


16± 5 


e2 south 


(-3.9,-2.9) 


0.74±0.03 


3.1±0.6 


9± 6 


e2 south 


(-2.0,-6.3) 


0.37±0.03 


6.5±1.4 


19± 6 


el 


(-4.2,-6.3) 


1.79±0.03 


2.1±0.3 


29± 4 


el 1mm peak 


(-4.2,-7.5) 


1.45±0.03 


3.0±0.3 


28± 3 


Peak of polarized flux 


(-6.5,-7.5) 


0.30±0.03 


5.8±1.6 


7± 8 


el 


(-2.0,-9.5) 


0.20±0.03 


10.0±2.7 


26± 7 


cl 


(-4.2,-9.5) 


0.63±0.03 


4.8±0.8 


28± 4 


cl 


cl average 


- 


3.2±0.6 


23± 5 




e2 average 


- 


1.1±0.3 


15± 7 b 




total average 


- 


1.8±0.4 


21± 6 b 





"Offsets are measured with respect to the phase center: Q2000=19' l 23 m 44?2, 52000= 
b e2 pol NW is excluded. 



14°30'33'.'4. 



Table 2 
Estimates of magnetic field strengths in the plane of the sky 



Region 


MT 


Hobs{°) b 


W (°) c 


riH ( cm 


3\d 


B p (mG)° 


el 


4.9±1.4 


4.8±2.0 


<5.8 


5xl0 4 
10 6 




> 0.3 

> 1.3 


e2 


6.7±3.4 


6.6±3.0 


<9.0 


5xl0 4 
10 6 




> 0.2 

> 0.8 



"The measurement uncertainty of the polarization angle 
b The observed polarization angle dispersion 
The intrinsic polarization angle dispersion 
d The number density of atomic hydrogen 
e The plane-of-sky magnetic field strength 
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Fig. 1. — Polarization map of the W51 el/e2 cores. The contours represent Stokes / at —5, —3, 3, 5, 10, 20, 30, ..., 110 a levels. The 1 a 
noise level of Stokes / is 27 mjybeam -1 . The line segments are polarization vectors, and their lengths are proportional to the polarization 
percentage with a scale of 3% per arcsec length. The dashed line that crosses the saddle point of the continuum is used to separate the data 
of W51 el and W51 e2 for our statistics. The star symbols mark the positions of e2, e4, e8, el and e3 from north to south. The triangle 
marks the position of H2O masers near e2 (W51 Main). e2 pol NW refers to the two vectors northwest of W51 Main. The ellipse indicates 
the synthesized beam of the Stokes Q and U images, which is 3'/2x2'/3 with PA=2°. Positions with high intensity but without polarization 
detections may have higher fraction of free-free emission or/and have an environment inappropriate for the magnetic alignment of the dust 
grains. 
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Fig. 2. — Distribution of polarization angle in W51 el/c2. The average and the dispersion of the polarization angle in el and e2 are labeled, 
and the double arrow indicates that e2 pol NW with <j> = 50° — 60° is excluded in the calculation. 
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Fig. 3. — The variation of the polarization percentage with Stokes / in W51 el/e2. The dashed lines are the results of the least-squares 
linear fit. 



